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^ Within the mSUGRA model, the observabihty of the decay of the next to hghtest neutrahno 

00 -^2 into leptons and the lightest neutralino Xi ^bls been studied using a full simulation of the 

CMS detector. The final state signature consists of two opposite sign leptons, several hard jets 
and missing transverse energy. The expected precision of the measurement of the dilepton 
CO mass edge is reported for 1 fb~^ of data, including systematic and statistic uncertainties, 

comparing two benchmark points with different signatures. 

o 

"O 1 Introduction 

> 

The standard model of particle physics (SM) provides no solution for pressing questions arising 
from astrophysical observations such as dark matter. In Supersymmetry (SUSY) a natural 
^ candidate for dark matter can be found if R-parity conservation is assumed. Supersymmetric 

particles (sparticles) have not been observed up to now what implies that they have to be heavy. 
On the other hand to provide a solution for the hierarchy problem their masses have to be in 
the TeV range. 

The long anticipated start of the Large Hadron Collider (LHC) in 2009 will allow to explore 
this new TeV range. With its center of mass energy of 10 TeV in 2010 it will allow to probe 
supersymmetric models very early on. A key point after discovery will be the determination of 
the sparticle properties. Because the lightest neutralino escapes detection, no mass peaks can 
be observed in SUSY decay chains. Of special interest are robust signatures such as mass edges 
in leptonic final states which can be probed with the CMS experiment'^. 



2 Leptonic decay of the next to lightest neutralino 



The leptonic decay of the next to lightest neutralino provides a clear signature of two opposite 
sign same flavour leptons which allows to trigger and identify the events in the hadron collider 



Table 1: mSUGRA benchmark points LMl and LM9. 
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Figure 1: Invariant mass distribution of leptons from the decay of the next to lightest neutralino are shown at 
generator level. The 3-body decay occurs at LM9 (a) and the 2-body decay at LMl (b). 



environment at LHC. In the framework of minimal Supergravity (mSUGRA) two leptonic decay 
modes of the next to lightest neutralino are possible a 2-body or 3-body decay. Two benchmark 
points (called LMl and LM9) have been studied which correspond to these decay modes. The 
mSUGRA parameters and the endpoint in the dilepton invariant mass of both points are listed 
in Tab. [TJ The low energy mass spectra of the two benchmark points have been calculated using 
the Softsusy code'^. Depending on the mass spectrum different decay modes are possible. 

A mass difference of the neutralinos smaller than the Z-boson mass and any slepton mass 
leads to a three body decay. In that case the endpoint represents directly the mass difference of 
the two lightest neutralinos 

mi, max = m^o - m^o. (1) 

The shape of the distribution depends on the mSUGRA parameters and is shown in Fig. [l] (a) 
for the LM9 benchmark point. 

A two body decay occurs via a real slepton and is allowed if one slepton is lighter than the 
mass difference of the neutralinos. In that case the endpoint can be expressed by 

(m~ — m?o ) [mlo — m?n ) 

) - ZI2 l^i 

m- 

The shape of the mass edge results only from kinematics and is triangular as shown in Fig. [T] (b) . 



3 Event selection 



In order to select the signal events a single leptonic trigger has been used. Because of the long 
cascade decays the final state consist of a high number of hard jets. The escaping neutralino 
leads to missing transverse energy in the detector. The selection requires three jets with > 
120 GeV, £:f > 100 GeV and > 80 GeV which are corrected for response. Additionally 
a missing transverse energy measured in the calorimeter (corrected for muons and jet energy) 





SUSY LMS 




tt+jets 




Z+jets 




W-f-jets 




{WW,WZ,ZZ)+jets 




QCD 




Background estimation 



CMS preliminary 




(3 

m 40 



10 



"(b) 






SUSY LM9 




tt+jets 




Z+jets 




W+jets 




(WW,WZ,ZZ)+jets 




QCD 



CMS preliminary 



50 100 150 200 2 



250 300 
m„ [GeV] 



Figure 2: The invariant mass distribution for same flavour opposite sign electron and muon pairs is shown (a). The 
opposite flavour opposite sign flavour background (b) consists mainly of top pair events and is used to measure 
the flavour symmetric background. The points in (a) represent the background extrapolation from (b). 



of at least 125 GeV is required. Finally two well identified and isolated {J2pt < 1-5 GeV for 
tracks within a cone of AR < 0.3) leptons with pip > 10 GeV are required. 

As main Standard Model backgrounds events with real isolated leptons are considered. To 
simulate these events the full CMS detector simulation has been used. The processes tt+jets, 
Z+jets and VK+jets have been simulated using the Alpgen MC-. The QCD background, which 
does not contain isolated high pT leptons has been simulated with the Pythia cod^. The SUSY 
sample is simulated using the Softsusy, SUSYHiflSI and Pythia codes. 

If one applies the full selection on the simulated set of events the main remaining background 
consists of tt+jets events Fig. [2] (a). Every background where the two leptons are produced 
uncorrelated, e.g. tt events, can be measured directly from data by selection of opposite sign 
opposite flavour lepton pairs as shown in Fig. [2] (b) . Using these wrong pairings one can predict 
the background distribution in the opposite sign same flavour lepton distribution. The method 
relies only on the knowledge of the lepton reconstruction efficiency which can be measured from 
data using events with a Z-bosons decaying leptonically'^. 

4 Measurement of the mass edge 

A combined fit consisting of three parts has been used to measure the mass edge. In case of a 
3-body decay the signal model consists of a quadratic term convoluted with a gaussian 

rricut 

f{x) = ^^ / dyy^e^^. (3) 



As background model a Landau curve has been fitted to the flavour symmetric background 
distribution Fig. |2] (b). The Z-peak is modelled by a Breit-Wigner function convoluted with a 
gaussian. The estimator of the mass difference {rricut) has been calibrated using MC templates 
with different values of the neutralino mass difference. The final fit for LM9 is shown in Fig. |3](a) 
and yields a value of 

{63A±l.lstat.±0.8syst.) GeV. (4) 

Similarly the endpoint in case of a two body decay has been studied ^. In that case the 
theoretical shape is known and a triangle, which is not calibrated using MC templates, is used 
to fit the distribution. The final fit shown in Fig. |3] (b) for LMl yields a value of 

{78.0±0.6stat.±0Asyst.) GeV. (5) 




Figure 3: Fit of the invariant mass distributions at LM9 (a) and LMl (b). The dashed green line in (a) is the 
background landau function, the dark green curve a Breit-Wigner function (Z contribution) and dashed red line 
the SUSY model. The dashed green curve in (b) is the background model and the red dashed line the SUSY 
model (triangle). The points represent the monte carlo events and the blue line is the final fit. 



The systematic error is calculated taking into account a misalignment and miscalibration 
of the CMS detector as expected after 100 pb^^ of data. Additionally the uncertainties on the 
jet energy scale (10%), the lepton energy scale and the acceptance have been included as well 
as uncertainties due to the fit model. It is found that this measurement is still dominated by 
statistical uncertainties at an integrated luminosity of 1 fb^^. The difference in the shape of 
both distributions can be used to identify the type of decay mode on a goodness of fit basis. 

5 Conclusion 

The leptonic endpoint can be reconstructed within the first LHC data (1 fb~^) if a low mass 
SUSY scenario is realised in nature. At both studied benchmark points the expected endpoint 
can be reproduced. The total error on the measurement is expected to be less than 3%. 
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